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ABSTRACT 

The physical state of interstellar gas and dust is dependent on the processes which heat and cool 
this medium. To probe heating and cooling of the ISM over a large range of infrared surface bright- 
ness, on sub-kiloparsec scales, we employ line maps of [C n] 158 /im, [O i] 63 /urn, and [N n] 122 
(Ltm in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard Herschel. We 
matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aro- 
matic hydrocarbons (PAHs) . We confirm at small scales in these galaxies that the canonical measure 
of photoelectric heating efficiency, ([C n] + [O i])/TIR, decreases as the far-infrared color, vf v (70 
/zm)/i//„(100 /im), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emis- 
sion, ([C n] + [O i])/PAH, is a near constant over a wide range of far- infrared color, 0.5 < 
vfv(70 /im)/z//„(100 /im) < 0.95. In the warmest regions, where vf v (70 fj,m)/i/f v (100 /im) > 0.95, 
the ratio ([C n] + [O i])/PAH drops rapidly to 4%. We derived representative values of the local 
UV radiation density, Go, and the gas density, rijj, by comparing our observations to models of pho- 
todissociation regions. The ratio Gq/hh, derived from fine-structure lines, is found to correlate with 
the mean dust-weighted starlight intensity, (U) derived from models of the IR SED. Emission from 
regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small 
grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 /im PAH 
ratio in regions that exhibit a deficit in ([C n] + [O i])/PAH, suggesting that small grains are ionized 
in these environments. 

Subject headings: galaxies: individual (NGC 1097, NGC 4559) — galaxies: ISM — ISM: lines and 
bands 



1. INTRODUCTION 

The life of gas and dust in galaxies is largely cyclic in 
nature. Interstellar gas and dust collapse to form stars, 
stars in turn process the material via nuclear fusion and, 
at the end of their lives, return a portion of this processed 
gas back into the interstellar medium (ISM), thereby al- 
tering and enriching the gas to be used in the formation 
of the next generation of stars. One key to mediating this 
process is the heating and cooling cycle of interstellar gas 
reservoirs. The physical conditions of interstellar clouds 
determine if those clouds collapse and form additional 
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stars. If gas is unable to cool sufficiently, star formation 
is quenched. Thus, the efficiency of heating and cooling 
cycles in neutral gas will affect the evolution of the gas 
and dust in the ISM. 

Neutral gas is thought to be heated primarily 
through photoelectrons that are eject ed from dust 
grains by incident UV radiation (e.g.. IWatsonl [T972t 
iTielens fc Hollenbachl [19851 IWolfire et all 1199.1 120031) . 
One of the key contributors of photoelectrons is thought 
to be the nearly ubiquitous polycyclic aromatic hydro- 
carbons ( PAHs) that emit strongly in the mid-infrared 
(IR) (e.g., fBakes fc Tielendll994HWeingartner fc Draind 
120011) . Small particles such as PAHs are excited by 
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UV photons and ar e subsequently radiatively de-excited 
(|Draine et alJ 12007). The vibrational bands emitted by 
PAHs should thus indirectly trace the ejection of photo- 
electrons, which are thought to dominate the heating of 
the neutral ISM. Ejected photoelectrons may heat a wide 
range of environments from dense molecular gas (T ~ 10 
- 50 K, Tin ~ 10 3 - 10 6 cm -3 ) to diffuse warm neutral 
gas (T - 5000 K, n H ~ 0.6 cm" 3 ) (|Drainell20lil) . 

On the other side of the balance between heating and 
cooling are the far-IR fine-structure lines that cool the 
heated gas. As it is relatively easy to excite (AE/k ~ 92 
K), all but the coldest gas is cooled by the collisionally 
excited [C II] 158 /im line. The abundance of carbon and 
the low ionization potential (11.26 eV) results in the [C n] 
fine-structure line at 158 /im being one of the bright- 
est emission lines se en in normal star forming galaxies 
(jLuhman et alj|1998f >. This line emission escapes the gas 
and dust allowing the gas to cool and thus provides a 
measure of the gas heating rate, assuming that the ISM 
has reached a local equilibrium. In addition to [C n] 
emission, the IR portion of the spectrum also hosts other 
key fine-structure lines (e.g., [N n] 122 /im, [O i] 63 /im, 
[Si II] 35 /im, etc.). This confluence of fine-structure lines 
and the PAH complexes in the mid- to far-IR makes this 
spectral range especially effective at mapping the heating 
a nd cooling of the ISM. 

IMalhotra et al.l (|2001[ ) examined globally averaged 
emission from far-IR fine-structure lines in a sample of 
60 star forming galaxies as part of the ISO Key Project 
Survey of Normal Galaxies. They showed that there is 
a decrease of line emission relative to continuum emis- 
sion in the far-IR in galaxies with warmer far-IR col- 
ors. This trend has persisted as investigations have mea- 
sured e mission from C + in more than 180 unresolved 
galaxies (jBrauher et al.l 12008! ). Mid-IR studies have also 
shown that the luminosity of PAHs decreases relative 
to the luminosity of larger dust grains, rather than re- 
main constant with d ust temperature, w hen looking at 
warmer far-IR colors (|Helou et al.lll991[) . This decline 
of PAH emission relative to large grains has been ex- 
plaine d by invoking inten se radiation fields that ionize 
PAHs (jTielens et alJll999f ). Alternatively, the depressed 
ratio of PAH emission to total infrared surface brightness 
(TIR) could be due to sampling different mi xtures of gas 
phase s within the large (~70") ISO beam (|Helou et al.l 
120011) . Recently, the unprecedented sensitivity of Her- 
schel has enabled the d etection line deficits in all far-IR 
line fine structure lines (jGracia- Carpio et alll2Qlll ). This 
has been interpreted as an increased competition for UV 
photons in dusty galaxies, as traced by the ratio of the 
far-IR continuum and the mass of molecular gas observed 
in a galaxy. Thus, significant progress has been made in 
understanding the heating and cooling of the interstellar 
gas of nearby galaxies. However, these efforts have been 
hampered by limited spatial resolution, spectral cover- 
age, and sensitivity. 

With the Photodetector Array Camera & Spectrom- 
eter (PACS) onboard the Herschel Space Observatory 
we are for the first time able to target the dominant 
cooling lines of neutral and ionized interstellar gas over 
sub-kiloparsec scales in nearby galaxies. The Kingfish 
project (Key Insights on Nea rby Galaxies: a Far- In- 
frared Survey with Herschel, iKennicutt et al.l 1201 lh is 



an open-time Herschel key program investigating critical 
heating and cooling diagnostics of dust and gas within 
nearby galaxies. In particular, Kingfish targets fine- 
structure lines that are the principal coolants of pre- 
dominantly neutral gas ([C II] 158 /im and [O i] 63 /im) 
and strong diagnostic lines of ionized gas ([O III] 88 /im 
and [N n] 122 /im) in nearby galaxies t hat span a wide 
bread th of key physical parameters (jKennicutt et al.l 
120111 ). To complete the picture of heating and cooling 
of the ISM in nearby galaxies, we combine the far-IR ob- 
servations of the Kingfish program with mid-IR obser- 
vations available fr om the Spitzer Infrared Nearby Galax- 
ies Survey (Sings. IKennicutt et al.ll2003f ). The spectral 
coverage and spatial resolution of these combined obser- 
vations allow us to investigate the distribution of PAH 
features and far-IR line emission on scales previously re- 
stricted to the Milky Way and Magellanic Clouds. 

NGC 1097 and NGC 4559 were selected as Kingfish 
Science Demonstration Phase (SDP) targets. NGC 4559 
is a late spiral type galaxy (SABcd) at a distance of ap- 
proximately 8.45 Mpc. It has super-solar metallicitjQ 
(12+lo g(0/H) = 8.810) and s hows no signs of AGN ac- 
tivity (|Moustakas et alj|2010h . In contrast, NGC 1097 
is slightly more enriched than NGC 4559 (12+log(0/H) 
= 9.09) and is a Seyfert 1 galaxy that contains a 
bright central ri ng of star formation approximately 2 
kpc in diameter Qy loustaka s et al.l 120101 : IHummel et alJ 
Il987t iSandstrom et al.l I2010D . At a distance of 19.1 
Mpc (jWillick et allll997l ). the central PACS pointing of 
NGC 1097 is one of the brightest regions in the Kingfish 
sample, when ranked by average TIR surface brightness. 

These two systems clearly probe a wide range of en- 
vironments: intense, concentrated star formation in the 
ring in NGC 1097 to nominal, distributed star forma- 
tion in the disks of both galaxies. This range of environ- 
ments gives a broad spread in far-IR color space allowing 
a robust baseline over which we may study the heating 
efficiency. Indeed, the Kingfish program has been de- 
signed to sample the widest possible range of star form- 
ing environments found in the nearby Universe. While 
NGC 1097 and NGC 4559 represent only a small fraction 
of the entire Kingfish sample, they span a representa- 
tive range of surface brightness and radiation fields for an 
exploratory study of the heating and cooling mechanisms 
operating in the ISMs of nearby galaxies. 

In this paper we present a broad investigation and 
comparison of properties of the ISM in NGC 1097 and 
NGC 4559 by combining data from the Kingfish and 
Sings programs. For a detailed investigation of heating 
and cooling in the central ring and e nds of the bar in 
NGC 1097 see also lBeirao et alJ (|2012f ). In §2 we present 
the observations and discuss the data processing. As 
both chopped and un-chopped PACS spectral data were 
acquired, these data are compared to each other as well. 
To accurately trace total PAH emission despite varying 
spectral coverage, we derive a method to combine lim- 

1 C ompared to the solar oxygen abundance of lAsplund et al] 
(2009), 12+log(0/H) = 8.69±0.05. 

2 We have adopted for thi s work the abundances deriv ed using 
the theoretical calibration of Kobuhiickv & Kcwlcv (2004). Adop- 
tion of the empirical calibration of Pilyugin & Thuan (2005) would 
yield subsolar oxygen abundances for NGC 1097 and NGC 4559. 
This choice slightly affects the magnitude, but not the direction, 
of the difference in metallicity between NGC 4559 and NGC 1097. 
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Fig. 1. — Observational footprints for IR spectral observations of NGC 1097 (left) and NGC 4559 (right) overlaid on MIPS 24/xm images. 
The green box denotes the area which had full Long-Low coverage with IRS. The magenta regions mark the area observed in Short-Low 
with Spitzer-IRS. The PACS-WS observational footprint is shown in Blue. PACS-CN coverage is shown in cyan. 



ited spectral and photometric data sets in §3. In §4 we 
discuss the heating and cooling of gas in these galaxies. 
The far-IR line deficit is investigated across the full range 
of star-forming environments in these two galaxies using 
the two strongest coolants, the [C n] 158 /zm and [O i] 
63 /Ltm lines. We discuss how the cooling is related to 
observed PAH emission in the ISM and investigate sce- 
narios that have been proposed to explain the far-IR line 
deficit in §5. Finally, we summarize the results of this 
study in §6. 

2. OBSERVATIONS 

2.1. Acquisition and Processing of P ACS Data 

Using PACS onboard Herschel we carried out far-IR 
spectral observations of NGC 1097 and NGC 4559 as 
part of the Kingfish SDP program. We report a detailed 
summary of observations in Tables [T] and [2] The integral 
field unit (IFU) of PACS has a roughly 47"x47" field of 
view composed of 25 9'.'4 x 9'.'4 pixels (fPoglitsc h et al.1 
2010). In Figure [1] we show footprints of the spec- 
tral observations of NGC 1097 and NGC 4559. These 
are aligned with Spitzer Infrared Spectrograph (IRS, 
iHouck et all l200l Long-Low (LL) coverage. The long 
strips measure 240" x 60" and 290" x 60" for NGC 1097 
and NGC 4559 respectively. In NGC 4559, observations 
are aligned with the major axis; whereas in NGC 1097 
the LL strip is nearly orthogonal to the star forming bar. 
The position angle of observations was determined by 
constraints in the pointing of the Spitzer telescope. In 
this work we highlight resolved mapping observations ob- 
tained for three key emission lines ([C 11] 158/xm, [N 11] 
122^m, and [O I] 63/xm) in both Chop-Nod (CN) and 
Wavelength Switching (WS) modes (see Tables [1] and 

121). 

All PACS spectral observations were reduced using 
the Herschel Interactive Processing Environment (HIPE) 
version 8.0 (|Ottl [2010V Reductions applied the stan- 
dard spectral response functions, flat field corrections, 

3 This now deprecated mode was superseded by the unchopped 
mode of PACS spectroscopy. 



and flagged instrum ent artifacts and bad pixels (see 
iPoglitsch et al.l l2010f ) . A dark was applied only to ob- 
servations taken in WS mode, as CN observations re- 
move the dark current via chopping. As inflight flux cal- 
ibrations are not yet available for PACS spectral data, 
ground based calibrations were applied to the data. 
These calibrations result in absolute uncertainties on the 
order of 30% and relative uncertainties of 10% and 20% 
in the blue and red grating orders respectively. 

Wavelength switching data were offered as an alter- 
native to chopping when a clean chop position was not 
available. In WS mode the grating is dithered by approx- 
imately half the FWHM of an unresolved line profile in 
an AABBBBAA pattern, where A is the reference grat- 
ing position and B is the grating position of interest. The 
flux at the grating position of interest may be compared 
to the flux detected at the reference position to remove 
the baseline, rather than using a reference at the same 
grating position with a spatial displacement. Observa- 
tions of an off position before and after each observa- 
tion supplement the spectral dithering. Following the 
standard reduc t ion pr ocedures for PACS, as outlined by 
IPoglitsch etaLl (|2010l) . we would subtract signal detected 
in the reference grating position from the grating posi- 
tion of interest. As the grating dither is small compared 
to the line width, this subtraction creates a differential 
line profile that must be modeled to reconstruct the data 
obtained in WS mode. 

Fortunately, because of the exceptional thermal sta- 
bility of PACS, we were able to process WS mode ob- 
servations without modeling the differential line profile. 
Instead, observations were processed under the assump- 
tion that the baseline was stable for the duration of the 
observation. This meant that WS reductions were not 
dependent upon models of line profiles to measure the 
total flux. However, due to drifts and long term tran- 
sients, the astrophysical continuum was undetected. In 
essence, this reduction process was the same as that used 
for the Unchopped mode that has since replaced WS in 
post SDP observations. 

The averages of the clean off observations obtained 
were subtracted from observations to correct for the ther- 
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mal background contributed by Herschel. Subsequently, 
all spectra within a given spatial element were combined 
after offsetting each to the median continuum value. Ap- 
proximately one fifth of the measured ramps exhibited 
a highly unphysical continuum resulting from transients 
in the detector. These transients were strongly corre- 
lated with motions of the grating and of Herschel it- 
self. Pixels which deviated from the mean spectrum by 
more than five standard deviations were masked out. Fi- 
nally, spectral cubes with 5" spatial pixels were created 
by projecting dithered raster positions of each spectral 
line together. In CN data, the two nod positions of the 
PACS array are slightly rotated on the sky with respect 
to one another. Therefore the nods are treated as sepa- 
rate raster positions when projecting the observations. 

To optimize the signal to noise and still take advantage 
of the spatial resolution, we employed 10" apertures tiled 
every 5" along the observed footprint. This aperture cor- 
responds to ~930 pc at the distance of NGC 1097 and 
~410 pc at the distance of NGC 4559. As apertures are 
not completely independent of their neighbors, there will 
be correlation between adjacent measurements. How- 
ever, this approach avoids introducing bias into the sam- 
ple by selecting apertures by eye. Furthermore, it fills 
out as large a parameter space as possible in the envi- 
ronment sampled for this study. 

Some apertures exhibited a residual baseline after pro- 
jection. To remove this, we fit a first order polyno- 
mial to the spectral range outside one FWHM on ei- 
ther side of the line center as estimate d from HI obser- 
vations of NGC 1097 and NGC 4559 (jKoribalski et al.l 
l2004HFischer fe Tullvl[l98lT) : this baseline fit to the spec- 
trum was then subtracted. Subsequently, line fluxes for 
each aperture were determined from PACS observations 
using two methods. We measured the flux by directly 
integrating under the continuum subtracted line. Addi- 
tionally, we used Gaussian profiles fit to the spectrum 
extracted from each aperture to determine the flux (see 
Figure [2]). For sources with detections above the three 
sigma level, we found that the flux measurements from 
both methods agreed within the relative uncertainties 
(20%) in NGC 4559. Lines from the central region of 
NGC 1097 are not well described by a Gaussian pro- 
file due to the large velocity dispersion of the star form- 
ing ring (see Figure [2]). Thus, the two flux determina- 
tions differed significantly in this spatial region. We have 
therefore adopted fluxes determined from direct integra- 
tion. Outside the ring, both methods of integration over 
a line are in agreement. 

The PACS spectral observations of NGC 1097 and 
NGC 4559 are shown as projected line maps in Figures 
[3] and |U Line maps of NGC 1097 have 3a detection 
thresholds at a surface brightness of ~ 0.04, 0.07, 0.01, 
and 0.07 x lO" 7 W m" 2 sr" 1 in [C n] CN, [C n] WS, 
[N n] CN, and [O i] CN lines respectively. Line maps 
of NGC 4559 have 3 c detection thresholds at a surface 
brightness of - 0.03, 0.05, 0.01, and 0.03 x 10~ 7 W 
m~ 2 sr- 1 in [C n] CN, [C n] WS, [N n] CN, and [O i] 
CN lines respectively. While the [N ii] 122 ^m and 
[O i] 63 /jm lines were only detected in roughly 50% of 
PACS pixels, even short integrations of the [C ii] 158 
/zm line yielded detections in 80% of the PACS pixels. 
Indeed, only when observations were taken off the axis 
of NGC 1097 was [C n] detection limited. 



2.2. Supplementary Data 

Both galaxies were al so observed by the Si NGS program 
(iKennicutt et al.l 120031 ) with Spitzer-IRS (|Houck et al.1 
|2004|). The central regions of NGC 1097 and NGC 4559 
were targeted by the short-low (SL) module. The two 
extra-nuclear regions targeted in NGC 1097 were also 
observed with SL. In addition to SL observations, the 
long-low (LL) module targeted strips across each target 
galaxy (see Figure Q}. There was no overlap of LL ob- 
servations with the extra-nuclear positions in NGC 1097. 
For a more detailed description of the Spitzer data, in- 
cluding reduction techniques and a presentation of the 
central s pectra of NGC 1097 and NGC 4559 we refer the 
reader to lSmith et al.l (|2007al ). For examples of the small 
aperture IRS s pectra and maps o f PAH emission we refer 
the reader to IBeirao et al.l (|2012f ) 

Spitzer-IRS spect ra were extracted using CUBISM 
( Smit h et al.l l2007bf) . Broad spectral features and the 
mid-IR lines were fi t and decomposed using PAHFIT 
([Smith et al.l l200"7al ). We fixed the silicate absorption 
optical depth at zero as ISmith et al. (2007a) found that 
negligible amounts of absorption from silicate dust were 
needed to reproduce the equivalent widths of lines in kpc 
scale regions of IRS spectral observations of NGC 1097 
and NGC 4559. Only data that had a surface bright- 
ness larger than three times the root mean square (RMS) 
noise were accepted as physically meaningful detections. 

The PACS and IRS spectral observations are supple- 
mented by imaging from 3-160 /im obtained as part of 
the Kingfish and Sings programs. These images are 
essential to the determination of the total infrared con- 
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Fig. 2. — Representative continuum subtracted spectra from re- 
gions with [N o] 122^tm detections in diffuse regions of NGC 1097 
(top) and the nuclear region of NGC 4559 (bottom). The [O i] 
63/jm line and the [N n] 122/mi line are plotted on the same scale. 
All spectra shown were extracted from the chop-nod observations. 
Dotted blue lines show Gaussian fits to each line. For spectra fro m 
the ring of NGC 1097 we refer the reader to IBeirao et al.l {2010). 
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Fig. 3. — Contour maps of far-IR fine-structure lines in NGC 1097 overlaid on the MIPS 24 (im image. The footprint of the spectral 
mapping in each line is shown as a faint purple outline. PACS beam sizes at the appropriate wavelengths are shown at the lower left of 
each panel. Contours begin at a SNR of 3<r; [C II] contours are plotted every 1.2 X 10~^W m -2 sr _1 , until 4 X 10 _8 W m -2 sr _1 ; [N II] 
contours are plotted every 1 X 10 _9 W m — 2 sr — 1 until 4; [O i] contours are plotted every 2 X 10 — 9 W m — 2 sr _1 until 1.4 X 10 _8 W m — 2 
sr . Line maps have been scaled to show the diffuse emission. For line maps of nuclear emission we refer the reader to[Bcirao ct ah 1 12010 ) 



tinuum emitted by dust, the far-IR color, and to aid in 
accounting for PAH emission where full spectral coverage 
is not available. The reduction of obse rvations followed 
standa rd procedures, as de tailed by iSandstrom et al.1 
IpOlOh and lDale et ail d2005h . 

2.3. Convolution of Data 

We convolve data to ensure that the comparison of 
fluxes data obtained at different resolutions is valid. This 
is particularly important as we are using a variety of 
instruments to cover a large wavelength range (3 - 160 
/im) . We are also interested in resolving the signatures of 
heating and cooling on small scales. These interests must 
be balanced as they lead us in diametrically opposite 
directions in the treatment of the data. 

Convolving to a larger point-spread function (PSF) will 
degrade the resolution of our data, reducing the dynamic 
range of heating properties being sampled. This is fur- 
ther complicated by the fact that data near the edge of a 
map cannot be properly convolved as there is a hard edge 
past which no data is obtained. Therefore, at the edges 
of a map the flux is not conserved. While map edges lie 
well beyond the regions of interest in our imaging data, 



this is not the case for the spectral maps. 

To test where we are limited in convolving the data, we 
have convolved masks of our IRS spectral maps (e.g., the 
value of each pixel recording data is set to one and pixels 
outside the footprint are set to zero) to the largest beam 
size used in this study (P ACS PSF a t 160 /im, FWHM = 
ll'/Ol) using sTinyTim (|Kristll2002D . We found this ef- 
fect to be modest when convolving Spitzcr-IRS LL spec- 
tra. Regions 22" wide at the center of LL strips (see 
Figure [TJ conserve approximately 98% of the flux, be- 
fore dropping to 70% of the appropriate flux at the edge. 
However, the Spitzer-IRS SL spectra data suffer a sub- 
stantial loss of flux; only 90% of the flux is conserved at 
the center of the spectral maps, dropping to 30- 5 0% at 
the edges. Furthermore, iPereira-Santaella et al.1 (|2010L 
see Appendix A) show that the PSF of Spitzer-IRS spec- 
tral maps are not well reproduced by modeled telescope 
PSFs, suggesting that the generation of maps has signif- 
icantly altered the PSF. 

As a compromise between matching resolution and re- 
taining reasonable spatial coverage without losing flux, 
we have convolved data with a beam smaller than the 
PACS 100 pan band (FWHM = 6'.'97) to this resolu- 
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Fig. 4. — Contour maps of far-IR fine-structure lines in NGC 4559 overlaid on the MIPS 24 /im image. The footprint of the spectral 
mapping in each line is shown as a faint purple outline. PACS beam sizes at the appropriate wavelengths are shown at the lower left of 
each panel. Contours begin at a SNR of 3<r; [C Ii] contours are plotted every 1 X 10 _s W m~ 2 sr — 1 ; [N Ii] contours are plotted every 4 X 
10 _10 W m -2 sr _1 ; [O i] contours are plotted every 4 X 10 -9 W m~ 2 sr — 1 . 



tion and bootstrapped data with a beam larger than 
the PACS 100 fim b and down to this res olution. Con- 
volution kernels from lAniano et al.l (|2011[) were used to 
convolve IRAC 3.6 fim (FWHM = 1"67), IRAC 8.0 /im 
(FWHM = 2"03), MIPS 24 fim (FWHM = 6'.'52), and 
PACS 70 ^m (FWHM = 5"41) images to the resolu- 
tion of the PACS 100 ^m band (FWHM = 6"97). As 
the delivered PSFs of Spitzer-IRS are extremely compli- 
cated due to the mapping str ategy employed by SINGS 
(|Pereira -Sant aeITa et al.lr2010f) . we have adopted a single 
Gaussian kernel appropriate for a central wavelength in 
each IRS module. Convolution of the [O i] 63/im data 
(FWHM ~ 5.3") was performed assuming a Gaussian 
PSF. 

To overcome the resolution limits for PACS 160 fim 
(FWHM = ll'.'Ol), we convolve PACS 70 fim observa- 
tions to the resolution of PACS 160 fim. A representa- 
tive 70 fim I 160 fim ratio was measured in 12" apertures 
centered on the associated sp e ctral ap erture. Similar to 
the approach in iSmith et al.l (|2007a[ ) , we assume that 
this color ratio is the same in both the 12" and 10" aper- 
tures. This ratio is combined with the photometry from 
the PACS 70 fim observations that are convolved to the 



resolution of the PACS 100 fim band to derive a surface 
brightness for PACS 160 fim at this resolut ion. We cal- 
culate the TIR following the prescription of iDraine fc Lil 

Hoo3), 

TIR = 0.95(vS u ) 8 + 1.15(vS u ) M + (uS u ) 70 + {vS v ) xm 

(1) 

using the matched 24 fim and 70 fim photometry com- 
bined with the bootstrapped PACS 160 fim, where the S^ 
are the surface brightness values in the indicated band- 
passes. While the [C n] 158 fim line falls in the PACS 
160 ^m bandpass, contamination from this line is less 
than 2%, therefore we do not correct for its inclusion in 
the calculation of TIR. 

Following our approach with the PACS 160 fim band 
described above, we have bootstrapped the [C ii] data 
to the resolution of the PACS 100/im band. Here we 
assume that the ratio of [C n]/^/„(70) does not change 
between 12" apertures, which sample the full 160 fim 
beam, and the 10" apertures. This ratio was combined 
with PACS 70 fim observations that were convolved to 
the resolution of the PACS 100 fim band to derive a [C Ii] 
surface brightness in the 10" apertures. Bootstrapping 
[C ii]/TLR and [C n]/PAH in this manner alters their 



Resolving the FIR Line Deficit 



7 




Fig. 5. — Flux measurements from data acquired in WS mode 
are plotted against the flux measurements from CN data. Line 
fluxes were measured from spectra extracted from 10" square WCS 
matched apertures. Measurements of the [C II] 158/im line are 
plotted as blue triangles and green squares, for NGC 1097 and 
NGC 4559 respectively; measurements of the [N II] 122/im line are 
plotted as magenta triangles and red squares, for NGC 1097 and 
NGC 4559 respectively. The solid line is a one-to-one correlation; 
dotted lines are indicative of the 20% relative uncertainty of PACS 
due to calibration uncertainty. 

surface brightness by <5%. The most significant changes 
occur in the warmest regions of NGC 4559, which lie at 
the south-west edge of the [C II] map next to a bright 
point source. The increased coverage of the 12" aperture 
is not completely contained by the [C n] map, thereby 
skewing the measured ratio of [C ll]/i//„(70). 

To ensure that the conclusions of this study are not af- 
fected by convolving to data to the PACS 100 /zm band 
rather than the 160 /im band, we have also convolved all 
our data to the resolution of the PACS 160 /im band. 
Using the coarser resolution and 12" apertures to ensure 
the PACS 160/zm beam is fully sampled, we reproduce 
the trends presented in this paper; however, the coarser 
resolution and larger apertures greatly reduce the num- 
ber of reliable apertures. 

2.4. Comparison of CN and WS Data 

Standard spectral observations with the PACS spec- 
trometer require chopping between the target and a 
blank space of sky at a high frequency to remove the 
thermal signatures of the telescope and to correct for 
transients. The PACS chopper is limited to a maximum 
throw of 6'. When targeting many nearby galaxies, an 
offset of 6' is insufficient to avoid emission from the target 
galaxy in the offset field. NGC 1097 and NGC 4559 have 
minor axis diameters of 6'. 3 and 4'. 4 respectively. These 
sizes provide clean off positions within the 6' throw of the 
telescope chopper, given the orientation of the observa- 
tions. Therefore, these galaxies provided an opportunity 
to verify that un-chopped data can reproduce the line 
flux measurements from standard chop-nod observations. 



Line data acquired in WS mode have ^40% less wave- 
length coverage than CN data due to the decreased num- 
ber of grating steps taken in the WS observations. The 
increased spectral resolution of the blue channel results 
in insufficient wavelength coverage to ascertain an accu- 
rate placement of the continuum. Therefore, we restrict 
our comparison of these two modes, CN and WS, to data 
acquired from the red photoconductor array, which has 
lower spectral resolution. Spectral windows were cen- 
tered on the nominal velocity of each galaxy. At the 
center of NGC 1097 the line-of-sight velocity differs by 
^500 km s -1 on the northern and southern edges of the 
rotating inclined star forming ring. This leads to wide 
double peaked lines in many of the PACS spectral pix- 
els. This rotational broadening of the line results in an 
increased uncertainty in the bright [C n] line as even 
less continuum is available for accurate placement of the 
baseline. 

[C n] 158/im and [N n] 122/im fluxes measured from 
CN observations are plotted against matched observa- 
tions in Figure [5] CN and WS data clearly agree within 
the ^20% relative uncertainty of the flux calibrations. 
Observations of the [N n] 122 fim line in NGC 1097 show 
the largest deviations from unity. This is most certainly 
due to the combination of the intrinsic faintness of the 
[N n] line and the uncertainty of the continuum place- 
ment in WS mode in the broad lines at the center of 
NGC 1097. 

While both CN and WS modes prove to be sensitive 
to detections of emission lines, we adopt CN data when 
available due to the (1) increased wavelength coverage 
that enables a robust placement of the continuum and 
(2) the increased depth of the CN observations of [N ii] 
and [O i] which, in many cases, had integration times 
that were larger than the comparable WS observation 
by a factor of two or more. As both observation modes 
provide suitable line detections, we make use of the [C n] 
WS data to extend the spatial coverage of both galaxies. 

3. TRACING PAH EMISSION 

Fully sampling the total PAH emission requires cover- 
age by both the LL and SL modules of the IRS. How- 
ever, a significant portion of our PACS coverage area 
is observed by only a single IRS module: LL 15-40 /im 
in the strips, SL 5-15 /im in the extra-nuclear regions. 
To expand our measurement of small grain heating effi- 
ciency to all regions with detections of [C n] and [O i] it 
is necessary to deduce the total PAH power via available 
data. 

To estimate the total PAH power with limited spec- 
tral coverage, we use the limited areas in which full 5- 
38 /Ltm coverage of PAH emission exists in NGC 1097 
and NGC 4559 to formulate measures of the strength 
of aromatic bands for apertures that have either LL or 
SL coverage, but not both. We define the total PAH 
power a s the combined inte nsity of dust features tabu- 
lated bv lSmith etafl (|2007al see Table 3). These features 
were measured in the complete combined IRS spectra. 
We selected major PAH emission complexes in LL and 
SL that are representative of the overall PAH emission 
to trace the total PAH power in regions where observa- 
tions of only one module were carried out. We chose 
the 17 /im complex in the LL spectra to represent to- 
tal PAH emission as this is the strongest feature covered 
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by the LL module (14-38 //m) (e.g., iSmith et al.|[200l 
I Werner et al.1 12004') . We select the 11.3 /zm complex as 
representative of PAH emission in the SL module (5-14 
fim). The 11.3 /im complex is neither the only nor the 
strongest PAH feature observed by SL. We selected this 
feature as it was a clean band that was subject to less 
uncertainty in the placement of the continuum, given the 
relatively weak silicate absorption features in these ob- 
jects. This is a requisite attribute of a representative 
PAH feature in SL as the uncertainty in the placement 
of the continuum underlying these bands is increased if 
LL data are not fit along with SL spectra. 

In addition to spectral measures of PAH, the strongest 
PAH band complex is included in the IRAC 8.0 /im band- 
pass. To increase the precision of using photometric data 
to determine PAH emission, IRAC 3.6 /jm data are used 
to su btract the stellar c ontribution to the 8.0/zm band 
(as in iHelou at "ail 120041 ). We correct the 8.0 /im IRAC 
band for the stellar component by scaling, and then sub- 
tracting, the 3.6 /im IRAC images, as determined from 
interpolating stellar SEDs for galaxies i n the Local Vol- 
ume Legacy survey (jMarble et a.l.ll2010h . 

PAH* Phot = [i/£„(8.Q) -0.24x^(3.6)]. (2) 

We note that this prescription leads PAHp^ to be 
~10% larger than the total power in PAH features. This 
may be due in part to the fact that we have not removed 
the warm dust continuum which is also present in the 8.0 
/tin IRAC images. Howeve r, it may also i ndicat e that the 
stellar correction given bvlMarble et all (120101) is on the 
low side (|Helou at aLll2004t iZibetti fc Grovesll201ll) . 

The coefficients of linear combinations of photometri- 
cally and spectroscopically measured PAH bands were 
adjusted to minimize the RMS scatter. As seen in Fig- 
ure |5] where we plot the ratio of the estimated PAH 
power (PAH*) and the true spectroscopic PAH power 
(PAHtot), the estimated total PAH* is found to be 

PAH* LL = 0.7472 x[PAH Phot + 3.25 x vS u (17)}, (3) 



and 
PAH 



SL 



0.497 x [PAH* Phot + 3.59 x ^(11.3)], (4) 



using either LL or SL data as well as IRAC data, where 
available. We note that PAH*^ and YAR* LL are unbi- 
ased tracers of the total PAH emission and are robust 
to changes in the small grain population, e.g., the ion- 
ization state of the PAHs or changes to the grain size 
distribution, unlike individual PAH bands. For instance, 
a clear trend with PAH ionization exists when the 11.3 
/im PAH feature alone is used to ascertain the total PAH 
power, see the middle panel in Figure |6] 

Regions near the edges of our IRS maps suffered a loss 
of flux in the convolution of Spitzer-IRS cubes. To re- 
tain these apertures in our sample, we measured both 
PAH Phot and the appropriate spectroscopic PAH, either 
PAH* LL , PAH* L , or PAH TO t- A fraction of PAH* phot 
equivalent to the fraction of lost flux, as measured from 
convolved masks of IRS modules, was added to the ap- 
propriate spectroscopic PAH measurement. This scaling 
resulted in a <5% correction in the total PAH power 
attributed to these regions. 

We find the 1 a scatter with respect to the mean is 
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Fig. 6. — We plot the ratio of five estimators of the total PAH 
power and the true PAH power measured from the full (5-38 /im) 
IRS spectrum as a function of the 7/11.3 (im PAH ratio in regions 
where full 5— 38^tm spectra are available, e.g., the nuclear regions 
of NGC 1097 and NGC 4559. 

IRS11.3, PAH2 L and PAH^, respectively. Thus, using 
both the spectral and photometric information together 
to measure the total power in PAH emission significantly 
reduces the scatter when only LL is available. While 
there is not a significant change in the scatter between 
IRS11.3 and PAH*j L , the combined data show much less 
dependence on the ionization state of the PAHs. We can 
further reduce the scatter in PAH* SL by including all the 
main PAH features present, however, this increases the 
correlation between PAH* SL and the ratio of the 7.7/11/3 
bands and is strongly dependent on the placement of 
the continuum, which is uncertain without LL observa- 
tions. 

4. HEATING THE ISM 

The efficiency of heating interstellar gas is controlled 
by the size and ionization state of the small dust grains 
that emit photoelectrons. The photoelectric efficiency is 
defined as the ratio of gas heating (via photoelectrons) to 
total dust heating (via ultraviolet and optical photons j3, 



gas 
dust 



(5) 



0.16, 0.21, 0.10, 0.15 and 0.09 dex for PAH 



Phot 



,IRSi7 



(Mochizuki 2004). Dust heating and gas cooling can be 
observed through emission in the infrared. Thermal ra- 
diation from radiatively-warmed dust peaks in the far- 
IR while the neutral gas, heated by electrons liberated 
from small dust grains, cools via collisionally excited fine- 
structure lines. The most prominent of the cooling lines 

4 This definition of the heating efficiency is based on observ- 
able properties and therefor e differs from the definition given by 
Ticlcns & Hollcnbach ( 1985), i.e., the ratio of energy that goes into 
gas heating compared to the far-ultraviolet radiation absorbed by 
small grains and PAHs. 
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are the [C I II 158 ^m, [O i] 63 /mi lines which domi- 
nate in low (iWolfire et all 119951 I2003D and high density 
(jKaufman et all 119991 12006) environments respectively. 
These two lines together constitute the bulk of the cool- 
ing in neutral gas. Indeed, the third most luminous neu- 
tral gas line, [Si II] 35 /im, contributes roughly an order 
of magnitude less to the cooling r ate than both [C ill 15 8 
/im and [O i] 63 /tm individually (jKaufman et al.l l2006). 

Observationally, it is assumed that the far-IR line-to- 
continuum ratio traces e PE such that 



J [C n] + 7 [Q i] 
TIR 



fa[C ii] + r ?[0 i]) £ ^' 



(6) 



where rj is the fractional contribution of each line to the 
total gas cooling, and (v^q n] + ^[O i]) — This im- 
plicitly assumes contributions from other cooling sources 
(e.g., [C i], [Si ii], CO, etc .) are negligible. Under 
this premise, iMalhotra et al.l ()200lD demonstrated that 
the global value of e PE decreases in galaxies that exhibit 
a warmer far-IR color, F„(60 /tm)/F„(100 /im). This 
drop, or deficit, in e PE can be attributed to a decrease 
in the rate of photoelectric heating, r gas , as dust tem- 
perature and star formation intensity increase. As the 
local UV radiation density, Go, measured in units of the 
local average interstellar field, increases relative to the 
particle density, uh, the effects of grain charging gain 
in importance. Therefore a common explanation for the 
decre ase of e PE has been t he ionization of small grains 
(e.g., IMalhotra et aljfeOOlD . Photoelectric heating be- 
comes less efficient as grains are ionized, resulting in a 
decrease in cooling line emission relative to the warm 
dust continuum. 

However, some studies have indicated that a high Go 
alone cannot produce the range of l ine ratios obse r ved in 
galaxies that exhibit a line deficit. iNegishi et al.l (|2001l ) 
found Go/riH to be independent of far-IR color, indicat- 
ing the mechanism producing the line deficit is related 
to the gas density. Alternatively, an excess amount of 
dust in the ionized gas could add to the far-IR contin- 
uum w hile suppressing [C i ll emission ; leading to the 
deficit (jLuhman et alJl2o5l I Abel et~all r2009). Indeed, 
recent results from Herschel showed that all far-IR fine- 
structure lines, independent of whether the line origi- 
nates in ionized or neutral gas, exhib ited a deficit relative 
to FI R at high values of Lfir/Mh 2 (jGracia-Carpio et"al"1 
120111 ): these results were apparent in a variety of galaxy 
types (e.g., LINERs, Seyferts, H II galaxies, and high-z 
galaxies). Yet it is not clear that the line deficit is inde- 
pendent of galaxy type. Still, other studies indicate that 
variations in PAH abundance can have significant effects 
on the the efficiency of heating and cooling in the ISM 
(jHabart et al.ll200l . Intense radiation fields are found 
in regions where PAHs are destroyed (jRubin et al.ll2009h . 
thus the existence of a line deficit in regions with a warm 
far-IR color may be tracing changes in the PAH abun- 
dances. 

Although great strides have been made in under- 
standing heating and cooling in nearby galaxies with 
the best techniques and spatial res olution available 
(|Brauher et al.|[200^lRubin et al.ll2009h . the use of global 
values can make it difficult to ensure that emission av- 
eraged together arises from similar physical environ- 
ments. iContursi et al.l (|2002t) noted that the integrated 



far-IR emission measured in global studies appears to 
be dictated by relatively few active reg ions that domi- 
nate t he signal within the beam, whereas M izutani. et al.l 
(2004) showed that overlapping photodissociation re- 
gions (PDRs) within an aperture lead to discrepant far- 
IR diagnostics. It has also been shown that the ob- 
served [C ii] line deficit is increased if emission that 
does not originate from the dense gas a round massive 
stars (e.g., non-PDR emission) is present (jLuhman et al.l 
2003). Given our sub-kiloparsec spatial resolution, we 
can confidently delineate regions dominated by diffuse 
and dense phases, although there will inevitably be some 
mixture of phases even at this resolution. 

In Figure [7] we plot the ratio of strong cooling lines to 
TIB0 as a function of PACS 70/100 /im emission; his- 
tograms on the right of Figure [7J show the distribution 
of e PE among regions cooler (black) and warmer (red) 
than z// iy (70)/^/„(100)=0.95 (see §5). We use the ra- 
tio vf v {7Q)/vfv (100) to trace the dust temperature as it 
is observationally easy to determine and enables a more 
direct comparison to previous studies of the far-IR line 
deficit. We note that i//„(70)/^/„(100) cannot be read- 
ily converted into Go due to contribution s of stochas- 
ticall y heated grains to the 70 /im flux (jlngalls et al.l 
120011) . While a robust measure of the radiation field 
requires d etailed fitting of th e spectral energy distribu- 
tion (e.g- IDraine fc Lill2007D . we find vf v (70)/vf„{lQQ) 
is well correlated with the star light intensity heatin g ob- 
tained from SEP models, see lAniano et al.l ([2012D and 
IDale et al.l (l20ll) . 

Over a significant portion of both NGC 1097 and 
NGC 4559, we find that emission from the two 
most important cooling lines, [C 11] 158 /im and 
[O 1] 63 /im, respectively account for an average of 
0.50±0.01%, 0.13±0.05% of TIR in cool dust regions and 
0.22±0.01% and 0.07±0.04% of TIR in warm dust envi- 
ronments. Si milar to previous studies of global measure- 
ments (e.g.. iLuhman et all 119981: IMalhotra et~aTl 120011: 
iGracia-Carpio et al. 1201 1|) . we find that the e PE of re- 
solved regions of the ISM in NGC 1097 and NGC 4559 
drops in areas characterized by warmer IR colors. This 
agreement is highlighted by the dash-dot lines shown in 
the top panel of Figure [7J that delineate the paramete r 
space filled by nearby galaxies in IMalhotra et all {2001). 

Even with the inclusion of [O 1] , the dominant coolant 
for dense environments, we find that the e PE is decreased 
in environments with more intense heating (see the bot- 
tom panel in Figure [7]). It is obvious that e PE is not 
uniform throughout a galaxy. Rather, local conditions 
substantially alter the processe s that heat and coo l the 
ISM. This is in agreement with iRubin et al.l (|2009l ) who 
found that, compared to diffuse gas, regions of intense 
star formation in the LMC (e.g., 30 Doradus) exhibit a 
s uppressed e PE . 

iBrauher et al.1 (|2008|) showed that e PE tends to decrease 
as a function of the far-IR color. This decline of line emis- 
sion relative to TIR is readily evident in regions within 
NGC 1097 whereas NGC 4559 is remarkably uniform. 
While NGC 4559 does not exhibit a strong decreasing 

5 Most investigations of the far-IR line deficit have used FIR 
band (42-122/jm) as a proxy for the IR continuum rather than 
TIR (3-1100/im). Use of this band does not change these trends, 
as TIR/FIR ~ 2. 
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Fig. 7.— Top: [C n] 158 fim / TIR plotted as a function of the far infrared color, vf v {70)/vf v (100). Middle: [O i] 63 /Urn / TIR plotted 
as a function of the far infrared color, v f v (J0) / v f v (VXi) . Bottom: ([C n] 158 /im + [O i] 63 /im) / TIR plotted as a function of the far 
infrared color, j;/„(70)/^/„(100). Even when accounting for the oxygen emission from denser regions, the photoelectric heating efficiency 
decreases with increasingly intense radiation. Regions in NGC 1097 are plotted as blue triangles whereas regions from NGC 4559 are 
plotted as green squares. The large grey diamonds show the weighted mean when data are binned. Histograms show the distribution of 
heating intensity for regions with gas warmer (red) and cooler (blac k) than ^ , / 1/ (70)/ i // 1 /(1 00)=0.95. The dash-dot lines in the top panel 
represent 80% inclusion curves of the galaxies from the ISO study of Malhotra et al. (2001). We remind the user that there exists a 30% 
uncertainty in the absolute calibration of PACS data when comparing these curves to the data. The mean error for each data set are shown 
to the lower left. Filled points indicate regions wherein the full IRS spectrum has been obtained. 



trend with far-IR color, this is likely due to the paucity 
of regions with warm dust that are observed. We also 
note that an offset appears to exist between these galax- 
ies that increases the measured dispersion in the line-to- 
TIR ratio. As the e PE of a gas cloud depends on the 
physical conditions of the gas a nd its location, we can 
conclude that the trend seen in Mal hotra et al.l (pOOl 
and lBrauher et all (|2008j ) is not simply the result of mix- 
ing different phases of the ISM on a global scale, but is 
either the result of true differences between galaxies or 
the mixing of phases of the ISM on sub-kiloparsec scales. 

Theoretical models of the photoelectric effect show 
that PAHs and other small grains dominate the 
production of photoelectrons in neutral gas (e.g., 



iBakes fc Tielensl [1991 iWeingartner fc Draind l200lh . 
This prediction has been supported by the strong correla- 
tion of [C il] emission and mid-IR flux in nearby galaxies. 
In addition to providing photoelectrons, PAHs simulta- 
neously emit photons in the mid-IR via vibrational and 
bending transitions. These mid-IR bands should provide 
a more precise view of heating and cooling than inte- 
grated thermal dust emission, which includes emission 
from large grains that are inefficient at heating the gas 
in the ISM. As PAHs are more directly coupled to heat- 
ing efficiency than TIR, they may offer a more robust 
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Fig. 8.— Top: [C II] 158 fim / PAH plotted as a function of the far infrared color, uf v (70)/uf u (100). Middle: [O i] 63 fim / PAH plotted 
as a function of the far infrared color, i//„(70)/v/„(100). Bottom: ([C II] 158 (im + [O i] 63 fim) / PAH plotted as a function of the far 
infrared color, v /„ (70) / v f v (100) , where PAH is traced by either PAH, PAHg^ or PAfl^ L , as available. Regions in NGC 1097 are plotted 
as blue triangles whereas regions from NGC 4559 are plotted as green squares. Filled points indicate regions wherein PAH was measured 
in the full IRS spectrum, while open points denote estimates of PAH from Equations 3 and 4. The large grey diamonds show the weighted 
mean when data are binned. Histograms show the distribution of heating intensity for regions with gas warmer (red) and cooler (black) 
than i'/ I /(70)/i//i,(100)=0.95. The mean error for each data set are shown to the lower left. 



and precise measure of the heating efficiency, e.g., 
J [C n] + J [0 i] . 

ip AH =( "[cm + "[oi]^ (7) 

where e' PE is the ratio of gas heating (via photoelectrons) 
to dust heating from PAHs, traced by emis s ion in the 
mid-IR. In the sample used bv lMalhotra et ail (|2001[ ). the 
globally averaged ratio of [C n] and PAH em ission shows 
rema rkably little scatter (a ~ 0.16 dex, IHelou et al.l 
120011 ). compared to the ratio of [C n] and the infrared 
continuum. Indeed, they found that PAH emission be- 
tween 5 and 10 /jm was directly proportional to the inten- 
sity of [C ii] line emission. If PAHs and emission from 
the [C n] fine-structure line are closely linked, heating 
by PAHs must decrease at high f„(60 /im) / f„(100 /im). 



Thus, IHelou et all (120011 ) posited that large grains be- 
come relatively more important to heating as the incident 
radiation field increases in intensity. 

5. LINE DEFICITS IN NGC 1097 AND NGC 4559 

We show the ratio of fine-structure cooling line emis- 
sion to PAH emission, as traced by either PAH, PAH* SL 
or PAR^ L , as available, in Figure [8] The dispersion in 
([C n] + [O i])/PAH (cr - 0.13 dex) is reduced com- 
pared to the ([C n] + [O i])/TIR (a ~ 0.20 dex), sug- 
gesting that ([C n] + [O i])/PAH is a more direct mea- 
sure of the heating efficiency in the ISM. When we focus 
specifically on the emission from small grains, neither 
NGC 1097 or NGC 4559 exhibit a gradual decline of 
([C n] + [O i])/PAH as a function of the far-IR color for 
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regions described by a far-IR color cooler than ~ 0.95. 
We also note that the apparent offset of e PE between 
NGC 1097 and NGC 4559 has disappeared. This lack of 
offset suggests that C + and PAHs are spatially correlated 
to a greater extent than C + and TIR or that the C + and 
PAH abundances are linked. This is supported by re- 
cent Herschel results which indicate that TIR is indee d 
influenced by old stellar populations (jBendo et al.ll20lil ): 
these populations are inefficient at heating the gas, but 
easily heat interstellar dust. 

Over much of the range in z//„(70 /im)/^/„(100 /xm), 
the line strength of both [C n] and [O i] are proportional 
to PAH emission, indicating the vibrational emission of 
PAH bands is directly proportional to the rate of pho- 
toelectron ejection from those same grains. However, in 
regions of the ISM characterized by warm dust, z//„(70 
/tm)/^/„(100 fj,m) > 0.95, we see a significant drop in 
both [C n]/PAH and [O i]/PAH, signaling that the cor- 
relation between the ejection of photoelectrons and vi- 
brational excitation of PAHs is not constant over the full 
range of bulk dust temperature (see histograms on the 
right of Figured]). A value of 0.95 in the ratio vf v (70 
//m)/i//„(100 /im) is representative of the color above 
which e PE drops and does not imply a single temperature 
that defines a break in the relation or physical conditions. 
As such, this color indicates where a bulk change in the 
physical conditions is seen; we have denoted this color by 
the vertical line in Figures [71 El and |9] 

We investigate, in these spatially resolved galaxies, sce- 
narios that have been proposed to explain the far-IR line 
deficit in global studies to ascertain the cause of this drop 
in e PE at warm far-IR colors. Notably we investigate the 
likelihood that self absorption, soft radiation fields, den- 
sity effects on [C n] line emission, and further ionization 
of PAHs can explain this occurrence in NGC 1097 and 
NGC 4559. 

5.1. Self-absorbed optically thick [C II/ 

[C n] and [O i] luminosities could be diminished 
if the photons were absorbed by cool foreground gas 
(|Luhman et al.|[T998f) . Self absorption has been observed 
in bo th the [O i] 63 /an line in Arp 220 ((Fischer et al.1 
HM% and the [C n] 158 /im line in observations of 
warm molecular gas in the Milky Way (|Boreiko fc Betd 
|1997|) . However, these occurrences of self absorption were 
seen in highly obscured regions; Arp 220, a nearby ul- 
traluminous infrared galaxy (ULIRG), has been shown 
to have unusually high optical depth (riooum ~ 5, 
Rangwala et al.1 1201 If) : while the PDR observed by 
Boreiko fc Betzl (|1997l) was located on the far side of the 
molecular cloud, with a minimum possible optical depth 
of 0.6 in the [C n] line. 

Using the ratio of Ha/H/3 from Ke nnicutt et al.l (|2009l ) 
we find an average extinction, (A v ), of 1.40 and 0.50 for 
NGC 1097 and NGC 4559 respectively. Looking only 
at the 20"x20" circumnuclear regions of each galaxy 
yields slightly higher am ounts of extinction, 2.34 and 2.22 
(jMoustakas et al.1 12010). However, these values are all 
considerably lower than the extinction seen in Arp 220 
where self-absorbed [C n] has been detected. It may also 
be noted that the extinction in the nucleus of NGC 4559 
is only 0.14 magnitudes lower than the extinction seen in 
the ring of NGC 1097. Despite this similarity, regions in 
the nucleus of NGC 4559 are not described by the warm 



dust temperatures seen in regions with a line deficit. 

Cold foreground diffuse clouds with A v ~ 1.0 could 
produce an optical depth of ~1.5 in the [C n] line center. 
These clouds would have to contain cold C + at the same 
velocity, within a few km s _1 , as the emitting cloud. 
A correlation in both velocity and position is unlikely 
between unassociated clouds. A chance alignment is ex- 
pected to be more likely in edge-on systems, where more 
clouds could lie in the line-of-sight. Contrary to this ex- 
pectation, the line deficit is primarily observed in NGC 
1097, a roughly face-on galaxy, and not as robustly in 
NGC 4559 which is inclined 64° to our line of sight. This 
scenario becomes even less likely given that the position 
of these clouds must correlate with an increasing dust 
temperature. 

Even so, we can determine an upper limit on the 
amount of C + self- absorption present by using the sili- 
cate absorption features at 9.8 and 18 fim as a surrogate, 
since the mid-IR continuum absorbed by silicates origi- 
nates in clouds that are likely to produce self-absorption. 
Scaling the silicate opacity will yield only an upper limit 
to C + self absorption as silicate absorption is insensitive 
to the velocity structure of the intervening clouds. The 
de-convolved [C n] linewidth in most regions is narrower 
than the PACS instrumental profile at 157 /im (FWHM 
= 240 km s _1 ); thus, line widths in our data are domi- 
nated by the instrumental profile. Observations of [C n] 
emission in a galactic star-forming region and its asso- 
ciated PD R show widths of 5 -10 km s _1 using HIFI on 
Herschel 's (jDedes et al.ll2010h . which has a velocity reso- 
lution of 0.7 km s _1 . For narrow [C n] emission, i.e., 4 - 
10 km s , the [C n] line becomes optically thick at N(H) 
= 3.65 - 6.25xl0 21 cm" 2 (|Russell et al.lfl980h . assuming 
N(C+)/N(H) = 1.6 x 1Q- 4 . Using the 10 /zm silicate ab- 
sorption profile from iLi fc Drainel (|2001[ ). these column 
densities would produce a minimum silicate opacity, T9.7, 
of 0.135 - 0.270. We find marginal detections of silicate 
absorption in only 13 apertures where full Spitzcr-IRS 
coverage exists. These regions, which are located on the 
western edge of the central ring of NGC 1097, yield (T9.7} 
= 0.06, with a maximum of T9.7 = 0.142. 

While some self- absorption by C + is possible, the T9.7 
yielded by fitting Spitzer-IRS spectra is not consistent 
with large amounts of self-absorption exist. We conclude 
that the existence of the far-IR line deficit in these galax- 
ies is not attributable to self absorption, which should 
introduce, at most, a drop in ([C 11] + [O i])/PAH by a 
factor of two. Nonetheless, we cannot eliminate the pos- 
sibility that intervening molecular clouds may exist that 
increase the scatter we see in ([C 11] + [O i])/PAH and 
([C 11] + [O i])/TTR. 

5.2. Soft UV radiation Field 

[C 11] emission in the local ISM will decrease after O 
and B stars have died, as photoelectons are not ejected as 
frequ ently in the softer radiation field (e.g. JSpaans et all 
1994), whereas the dust, including PAHs, continues to 
be heated. The existence of a soft radiation field in the 
warmest parts of a galaxy seems unlikely, e.g., a high 
vfv(70)/vf„(100) does not imply an absence of O and B 
stars. Indeed, recent observations of Paschen a emission 
in the central ring have shown that the very regions that 
exhibit a line deficit have a higher star formation rate 
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Fig. 9.— Top: [N II] 122 /an / TIR plotted as a function of the 
far infrared color, i/f v (70)/i/f v (100). Bottom: [Si n] 35 /im / TIR 
plotted as a function of the far infrared color, uf v (70)/uf v (100). 
If a dusty H II region was the cause of the line deficit, then all lines 
that become optically thick should show the deficit, which is not 
seen in either [N II] / TIR or [Si II] / TIR. Colors and shapes of 
points follow the same scheme as Figure [8] 



than other parts of NGC 1097 dHsieh et al.ll2011l). The 
eleva ted star formation and FUV flux (Kuchinsk i et al.l 
120011 ) associated with the ring in NGC 1097, which is 
associated with many of the apertures that show a line 
deficit, lead us to reject a soft radiation field due to a 
dearth of massive stars. 

5.3. Dusty H II Regions 

The presence of large amounts of dust in the vicinity 
of photoioni zed gas could also lea d to the creation of the 
line deficit ([Luhman et al.l [2003) . Direct absorption of 
these photons by the dust also leads to an increase in TIR 
while decreasing the UV photons that escape and heat 
the neutral gas , thus leadi n g to a net reduction in e PE . 
The models of lAbel et"afl ([2009T I have shown that the 
presence of dust in regions with a high ionization param- 
eter removes Lyman continuum photons that would nor- 
mally heat the gas. The absence of these photons leads 
to a reduction of far-IR fine-structure lines from both the 
ionized and neutral gas. However, with inclusion of ra- 
diation pressure, large grains may be compressed into an 
ionized shell, leading to an increase in the stellar photons 
that ionize H, compared to a uniform dust distribution, 
negating the ability of dust to cause a decrease in far-IR 
line emission tiDraindl20lih . 

Neither [N n] 122 /jm /TIR nor [Si n] 35 /im / TIR 
show a trend with far-IR color, as illustrated in Figure 
[9] The lack of a deficit in the [N n] line, which is only 
found in ionized gas, is particularly difficult to explain 
by invoking a dusty H II region. If any trend is no- 
table in [N n]/TIR it is an increase in regions charac- 
terized by warm far-IR color. The different [N Ii]/TIR 
plateaus may be related to the metallicity difference 



between NGC 1097 and NGC 4559, as secondary pro- 
duction of nit rogen may exist at the higher metallicity 
of NGC 1097 ([Vila Costas fc Edmundslll99l . Alterna- 
tively, this effect could also be the result of H n regions 
in the ring of NGC 1097 filling a larger fraction of the 
ISM than other regions where [N n] emission is detected. 

Global studies of ULIRGS have indicated the pres- 
ence of dusty regions and a deficit in the [N il] 12 2 
fj,m line (jFischer et alj|2010t IGracia-Carpio et al.ll201lD . 
We see no indication of this effect in resolved regions 
of NGC 1097 and NGC 4559. While we cannot rule out 
the applicability of dusty H II Regions to ULIRGs, we do 
not find evidence of elevated dust lev els in the ionized gas 
within these two galaxies. However, IBeirao et al.l (|2012[ ) 
find L FIR /M H2 ~ 50 Lq/Mq in N GC 1097. While this 
is belo w the ratio of Lq/Mq where IGracia-Carpio et al.1 
see a line deficit begin (80 L^/Mp,), it lies near 
this d omain. If we plot data from IGracia-Carpio et al.1 
(|2011[) with the ring of NGC 1097, we find that regions 
on the central ring of NGC 1097 lie at the lower edge of 
the scatter in their data and transition nicely into the 
deficit seen in galaxies with higher Lfzr/Mh 2 - 

5.4. Contributions to [C II] from Non-PDR Gas 

To understand the far-IR line emission, we must also 
understand the conditions of the ISM that lead to far- 
IR line emission. We determine the average gas density, 
nn, and radiation intensity, Go within each region of 
NGC 1097 and NGC 4559, assuming that all the emis- 
sion originates from PDRs, by comparing [C n], [O i], 
and TI R in each aperture t o models of ionization from 
PDRs (|Kaufman et al.1120061 ) . Data are overlaid on plots 
of model [C n] / [O i] versus ([C n] + [01])/ TIR in 
Figure [U 

According to PDR models, the gas which dominates 
the cooling in NGC 1097 and NGC 4559 lies at densi- 
ties between 10 2 5 and 10 3 cm -3 . If the line deficit is 
caused by density effects, the [O i] line should become 
the dominant cooling line in the densest clouds. To the 
contrary, as shown in Figure [S] rather than becoming the 
dominant line, [O i] exhibits a deficit congruent with the 
[C n] deficit. This indicates that the observed line deficit 
is not a product of increased gas density. 

The total integrated [C n] emission along a line of sight 
includes contributions from ionized gas in addition to the 
PDRs. This additi onal [C ill em i ssion must be removed, 
as the models of Kaufman et ail (|2006T ) only account for 
emission from neutral gas. The [C n] contribution from 
ionized gas can be traced by emission from N + , which 
in the Galaxy is emitted p rimarily in the diffuse ionized 
phase (jBennett et al.lll994f ). 

This correction is straightforward given measurements 
of [N n] 205 /im, as this faint line has the same criti- 
cal density as the [C n] 158 /im line in ionized gas (~ 
45 cm" 3 ) (jOberst et al.ll2006l ). This results in a [C n] 
158 /im/[N n] 205 /mi ratio that is nearly constant (see 
Figure Hip , varying by only a factor of 1.6 over three 
orders of mag nitude in density ([Oberst et alJ l2006). Un- 
fortunately, the faint [N n] line at 205 /im is difficult to 
detect. The only observations of this line for these two 
g alaxies cover t he sta r forming ring in NGC 1097. 

IBeirao et al.l (|2012| ) employ the region in which [N n] 
205 /im emission is observed to determine that emission 
from diffuse ionized carbon is equal to ~1.2 x [Nil] 122 
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([c II] + [o I]) / TIR 

Fig. 10.— [C n] 158 /xm/[0 i] 63 /jm plotted against ([C n] 
158 (im + [O i] 63 /mi) /TIR in NGC 1097 (blue triangles) and 
NGC 4559 (green squares). A grid of Go, in units of the local 
Galactic radiation field, and njj, in cm" 3 , as determined from 
the PDR models of Kaufman ct al. (2006) are over-plotted. Only 
data with a reliable [N n] 122 /im line detection are plotted when 
correcting for an ionized component. In the middle plot we assume 
a n e of 100 cm -3 when the [S in] ratio is in the low-density limit. 
In the bottom plot we assume a n e of 10 cm -3 when the [S III] ratio 
is in the low-density limit. Error bars do not include uncertainties 
in the correction for diffuse [C II] emission. 
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Fig. 11. — The ratio of [C o] / [N II] from diffuse ionized gas 
is plotted as a function of n e for both the 122 /mi and 205 (im 
transitions. The range of [C n] / [N II] 122 /tm we expect to see 
in NGC 1097 and NGC 4559, based on density variations, is high- 
lighted by t he cyan strip. Fo r comp arison, the single value correc- 
tion used by Malhotra ct al. (2001) is shown as a diamond on the 
curve tracing [C II] / [N II] 122 /mi; whereas the circle denotes the 
[C II] / [N II I 122 /tm determined from the [N n] 122 /xm/[N n] 205 
/im ratio bv lBeirao et al.l p012T l. 

/im. As we are interested in studying emission extended 
well beyond the brightest regions we must allow for this 
ratio to change. There are no observations of [N n] 205 
/im outside the nucleus of NGC 1097 that allow us to 
track these changes. Furthermore, no observations of 
[N n] 205 /im exist to verify that this ratio is the same 
in NGC 4559, whic h has a lower star formation rate 
([Calzetti et al.ll2010D 

An alternative approach is to use the relatively 
stronger [N II] 122 /im line as a means of tracing the 
diffuse ionized [C ill emission (e.g.. lKaufman et"al1ll999t 
iMalhotra et al.ll2001l and references therein) . The higher 
critical density of [N n] 122 /im (~ 300 cm -3 ) leads to a 
stronger dependence of [C II] 158 /xm/[N n] 122 /im on 
density compared to [CII] 158 urn/ [Nil] 205 um. This 
is shown in Figure HT1 which plots the [C n]/[N n] line 
ratios as a function of electron density. Preferably, this 
density would be calculated from the ratio of [N n] 122 
/im/[N n] 205 /im but as discussed, the weaker line is 
unobserved for most of this sample. A substitute density 
indicator is available in the mid-IR ratio of [S III] 18.71 
/im /[S ill] 33.48 /im. 

The [S III] ratio is sightly temperature dependent and 
is most sensitive to d ensities in the range 100 - 10 4 cm~ 3 
(Sni iders et al.l [20071) . We measure this ratio in our LL 
coverage of NGC 1097 and NGC 4559, when these lines 
are detected, and derive an elec tron density, n e , using 
the models of ISmith et all (120091) . We derive [C n] 158 
/zm/[N II] 122 /xm appropriate for ionized gas as a func- 
tion of n e , in order to scale the [N n] emission appropri- 
ately. To determine this ratio, we used observed Galactic 
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Fig. 12. — The fraction of [C n] emission attributed to PDRs is 
shown as a function of vf v (70) / vf v (100). All methods of correcting 
for emission from diffuse ionized gas show a decrease in the PDR 
fraction at warm colors. Rather than eliminating the line deficit, a 
plot of [C u]pdr I PAH would significantly increase the line deficit 
as the fraction of [C II] emanating from PDRs decreases in warm 
regions. 

gas p hase abundances of C (1.6xl0~ 4 per hydrogen nu- 
cleus, [SofiaetanJlOOlLand N (7.5xlCT 5 per hydrogen 
nucleusT TMever et al.lll997f ). F or nitrogen we adopte d the 
transition probabilities from | Galavis et al.l (119971 ) and 
collision strengths from Hudson fc Belli (120051); while for 
carbo n we use transition probabilities f romlGalavis et ah 
(1998) and collision strengths from iBlum fc Pra dhan 
(1992J). 

The correction for diffused ionized gas using the [S in] 
lines to trace the density is in agreement with the cor- 
rection that would be adopted using the measurements 
of [C n] /[N n] 205 fim in the nucleus of NGC 1097, see 
Figure [TT] We show the magnitude of this correction 
visually by plotting our data on models of PDR emis- 
sion both with and without correcting for diffuse ionized 
gas in Figure [TUl There is a significant shift in 20% - 
40% of the data points with a 3cr detection of [N ii] 122 
/im, depending on assumptions in the low density limit, 
when emission from ionized gas is removed. Given the 
uncertainties in correcting for contributions to [C ii] from 
ionized gas, reality most likely lies between these repre- 
sentations of the UV field and the gas density we show 
in Figure ITOl 

When the [S III] ratio is in the low-density limit, n e < 
100 cm" 3 , the ratio of [C n]/[N ii] 122^m in ionized gas 
can vary between 1.24, at n e of 100 cm" 3 , and 6.8, at 
n K of 0.1 cm~ 3 . O bservations of the [N ii] 205 /xm line 
(|Beirao et al.l 120121 ) indicate that at least some of the 
regions that are in the low density limit have an n e of 
~ 100 cm" 3 ; therefore, in the nucleus of NGC 1097 we 
adopt the conservative scale factor of 1.24 appropriate 
for n e = 100 cm" 3 . In regions outside the central star 



forming ring of NGC 1097, the flux of the [N n] 122 fan 
line is sufficiently weak relative to [C ii] that the effects 
of this uncertainty are negligible (e.g., <5%). 

If observed PAH emission only receives contributions 
from the diffuse ionized gas, then e' PE should be cal- 
culated using only [C ii] emission that originates from 
PDRs. We show the fraction of [C n] emission that 
emanates from PDRs as a function of the far-IR color 
in Figure [12J In addition to the correction adopted in 
this work, we also calculate the fraction of [C ii] that 
would ori ginate in PDRs using t he single value prescrip- 
tion from iMalhotra et~a l. (2001). It is noteworthy that 
all regions show a relatively high PDR fraction, regard- 
less of the far-IR color. Furthermore, the fraction of [C ii] 
emitted by the PDR decreases in the warmer regions, re- 
gardless of which correction for ionized gas is employed. 
An gp E associated with pure PDR emission would thus 
exhibit an increased line deficit at warm far-IR colors, 
accentuating the sudden decrease in ([C n] + [O i]) / 
PAH. 

In addition to an ionized component and PDRs, the 
ISM also contains diffu se HI. Observations of the dust 
SED in these galaxies (jAniano et al.l 120121 ) are incom- 
patible with the high Go we determine assuming that all 
the [C ii] emission and TIR originates in a PDR with an 
implied starlight intensity heating, (U), that is inconsis- 
tent, typically 10 - 50 x lower than expected from Gq. In 
reality, some fraction of both TIR and the [C ii] line will 
emerge from the diffuse ISM. Attributing 70% of TIR 
and 30% of [C n] surface brightness to cool diffuse gas 
would bring the Go and from PDR model and (U) from 
the dust SED into closer agreement, typically (U) 1 - 
10 x lower than expected from Go- This will simultane- 
ously increase the nn of the modeled PDR, as shown in 
Figure fl~3l Standard models of starlight heating indicate 
that this diffuse gas is too cool to emit [O i]. Decreas- 
ing the amount of [C n] relative to [O i] emitted by the 
PDR could result in [C ii] no longer being the dominant 
coolant. However, the inclusion of other line coolants, 
e.g., [O i], does not eliminate the line deficit. It is note- 
worthy that models of starlight heated diffuse gas also fail 
to produce t he warm H2 emissi on seen in some translu- 
cent clouds (Ingal ls et al.ll20lH ). H 2 emission lines are 
seen in our Spitzer IRS spectra, indicating that the dif- 
fuse cloud population may also harbor hot material that 
supplements the [O 1] emission, requiring non-standard 
heating processes that are not fully described by a PDR 
model. 

Assuming that diffuse gas is distributed in a relatively 
uniform manner throughout the ISM, correcting TIR, 
[C 11], and [O 1] for its contributions would alter the ab- 
solute values of n e and Go, but would not change the gen- 
eral trend of points relative to each other, see Figure 1131 
PDR models are not able to reproduce all the spectral 
line and dust SED features we have observed. However, 
none of the corrections we have investigated would alter 
the fact that the ratio of Go/n e and (U) both increase 
in regions where the line deficit is observed. Contribu- 
tions of line emission and TIR from diffuse gas, both 
ionized and cool, cannot produce the far-IR line deficit 
in NGC 1097 and NGC 4559. 
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Fig. 13. — Similar to the middle panel of Figure flOl except that 
we assume 70% of TIR and 30% of [C II] emission originates in 
diffuse clouds. 
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PAH 7.7 / 1 1 .3 

Fig. 14.— Top: [C n] /PAH is plotted as a function of the PAH 7.7 
/ 11.3 /^m ratio. An increase in this PAH ratio, which is observed 
in all regions showing a deficit, is indicative of an increase in PAH 
ionization. Data shown here are limited to where SL observation 
had been acquired, see Figured Colors and shapes of points follow 
the same scheme as Figure [8) Bottom: Histogram of the PAH 
7.7 / 11.3 fim ratio. Regions for which ^/„(70 fim) / v /„ (100 fim) 
< 0.95 are indicated with the black histogram, whereas regions 
with f/„(70 /im) / v /„ (100 (im) > 0.95 are indicated with the red 
histogram. 



5.5. Grain Charging 

Positively charged dust grains are inefficient at heating 
gas due to the greater potential w ell that must be over- 
come to eject additional electrons (|Tielens fc H ollenbach 
1985). This decrease in gas heating is not accompanied 
by a strong decrease in the thermal dust emission that 
dominates TIR; nor does it cause the vibratio nal PAH 
emission to drop in strength (jDraine fcTHl2007l) . There- 
fore, an increase in grain ionization will lead to a deficit in 
line c ooling relative to t he thermal dust and PAH emis- 
sion (jTielens et al.lll999fl . 

Figure [TOl shows that most of the gas in NGC 1097 and 
NGC 4559 lies at a density between 10 2 5 and 10 3 cm~ 3 . 
However, the local UV field, Go, varies considerably, 
ranging from 50 - 1000, when the most conservative 
corrections for ionized gas are employed. We see the 
largest variation of UV heating in the ring of NGC 1097, 
while NGC 4559 displays remarkable uniformity in both 
Uh and Go across its disk. We compare ([C n] + 
[O i])/PAH and ([C n] + [O i])/TIR to G /n H de- 
rived from models of PDRs and the mean dust-weighted 
starlight intensity, {U), deri ved from dust models of the 
IR SED (jDraine fc Lill2007t) . where the average value of 
(U ) in each regi o n wa s extracted from maps produced 
by lAniano et al I (|2012f ). Both gas and dust diagnostics 
clearly indicate that strong radiation fields are incident 
on the regions showing a line deficit. Indeed, Gq/uh 
increases more than an order of magnitude between re- 
gions of low and high incident radiation, suggesting that 
rain charging could be important in these warm regions 
Wolfire et al. l [l990l) . 



If the ionization of PAHs is indeed responsible for 
the observed line deficit, a significant variation in the 
7.7 fim to 1 1.3 fim PAH band flux ratio should also 
be observed (jDraine "&HI20071) . The 7.7 and 11.3 
PAH features are both observed using the SL mod- 
ule and have very similar PSFs, which are consider- 
ably smaller than the 12" apertures we have used. As 
the convolution of Spitzer-IR S maps are very uncertain 
(jPereira-Santaella et al.l 120101 ) we compare the 7.7 fim 
to 11.3 fim PAH band flux ratio measured from un- 
convolved spectra with the same ratio measured from 
convolved spectra. As expected, convolution of SL spec- 
tra decreases the spread seen in this PAH band ratio. 
The un-convolved spectra clearly show that more ion- 
ized PAHs are found in regions where the line-to-PAH 
deficit is strongest. 

In Figure [14] we plot convolved [C n]/PAH as a func- 
tion of the ratio of the un-convolved 7.7 fim/11.3 fim 
PAH bands. Regions that show a far-IR fine-structure 
line deficit relative to total PAH emission do indeed ex- 
hibit a modestly greater than average 7.7 fim / 11.3 fim 
ratio, suggesting that PAH ionization fraction is higher 
in regions with a line deficit. This is consistent with an 
elevated Gq/uh seen in these regions. 

The shift seen in the 7.7 fim to 11.3 fim PAH band 
flux ratio could also be indicative of a change in the 
grain size distribution. Small PAHs are more efficient at 
producing photoelectrons, compared to the larger PAH 
grains, as th ey tend to stay neutra l rather than becom- 
ing charged (Bak es fc Tielenslll994D . This has been con- 
firmed via Galactic observations of the Horsehead neb- 
ula that show a smaller PAH/VSG (very small grain) 
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ratio in the warm gas compared to the diffuse ISM 
(|Compiegne et al.l 12008). However, we do not detect a 
correlation between either the 7.7 jtm/17 jtm nor the 11.3 
jim/17 fim PAH band ratio and ([C n] + [O i])/PAH. 
If grain size was the driving factor in the suppression of 
([C n] + [O i])/PAH, we should detect a significant trend 
with both of these ratios as the 17 /im complex is at- 
tributed to the largest PAHs. Furthermore, we find that 
QC n] + [O i])/PAH does not correlate with the 6.2/7.7 
/xm band ratio, which has traditio nally been associated 
with grain size (jDraine fc Lill200ll ). 



6. CONCLUSIONS 

In this paper we have sought to understand the phys- 
ical conditions that drive the heating and cooling of gas 
and dust in the ISM of two nearby galaxies using spectral 
mapping data of NGC 1097 and NGC 4559, which are 
among the first galaxies observed by the PACS spectrom- 
eter on the Herschel Space Observatory as part of the 
Kingfish program. We matched these maps to Sings 
mid-IR observations to extract spatially resolved infor- 
mation on both gas and dust heating in these two nearby 
galaxies. Using these combined observations we have in- 
vestigated photoelectric heating efficiency and the impor- 
tance of PAHs as a source of photoelectrons over a large 
range in surface brightness with higher spatial resolution 
than previously possible in external galaxies. 

As seen in other galaxies, we find that the photoelec- 
tric heating efficiency, as traced by ([C il] + [O l])/TIR, 
is lower where the dust temperature is elevated. Con- 
trary to a decreasing heating efficiency, which is evident 
in NGC 1097, ([C n] + [O l])/PAH is flat as the temper- 
ature of the dust increases, up to 1^(70 jim)/f/„(100 
jim) ~ 0.95. In this regime, gas heating, as traced by fine- 
structure cooling line emission, and dust heating, dis- 
cerned by the power emitted in the mid-IR from PAHs, 
effectively trace each other, indicating that the ratio of 
([C il] + [O i])/PAH is uniform in many environments. 
The scatter observed in ([C n] + [O i])/PAH is also re- 
duced compared to ([C n] + [O i])/TIR. This suggests 
that PAH emission delineates the ejection of photoelec- 
trons more effectively than the integrated dust emission, 
TIR. In gas with dust warmer than vf v {70 jtm)/i//„(100 
jtm) ~ 0.95, however, we observe a deficit in [C ii]/PAH 
and [O l]/PAH, that was not previously detected. 

Many explanations have been proposed to explain the 
deficit in the observed heating efficiency. Comparing 
the properties of regions in which this deficit is seen in 
NGC 1097 and NGC 4559 we find that: 



• Both [C n] and [O i] exhibit line deficits relative to 
both TIR and PAH emission. However, deficits are 
not observed in the [N n] or [Si il] lines, as would 
be expected if dusty Hn regions were responsible. 



Abel, N. P., Dudley, C, Fischer, J., Satyapal, S., & van Hoof, 

P. A. M. 2009, ApJ, 701, 1147 
Aniano, G., Draine, B. T., Gordon, K. D., & Sandstrom, K. 2011, 

PASP, 123, 1218 
Aniano, G. et al. in prep. 



• If the line deficit was the result of increased density, 
cooling by [O i] should compensate for the dearth 
of [C n] emission in dense gas. This is not the case, 
as [O i] and [C il] exhibit similar deficits. 

• Correcting [C n] emission by removing contribu- 
tions to this line from both ionized and diffuse gas, 
does not remove the line deficit. To the contrary, 
this correction would accentuate the [C n]/PAH 
deficit observed in regions with warm far-IR col- 
ors. 

• Both (U) from SED models and the ratio of Gq/uh 
from models of PDRs indicate that regions where 
a deficit is seen have intense radiation fields that 
could lead to grain charging. 

• The distribution of the 7.7 jim/11.3 jim PAH band 
strength ratio differs between warm regions, that 
exhibit a line deficit, and cool regions that do not. 
This indicates that grain charging, which leads to 
inefficient photoelectric heat, could be responsible 
for the observed line deficits. 

While NGC 1097 and NGC 4559 have regions spanning 
a large range of IR surface brightness, they represent 
a small fraction of the diverse environments probed by 
Kingfish. An analysis of the heating and cooling in 
the full sample will probe more intense starbursts, larger 
ranges of metallicity, and different types of ionization 
sources. By analyzing the full Kingfish sample we will 
be able to more fully study the range of gas and dust 
heating explored here and expand our investigation to 
even more extreme environments to test the universality 
of these results. 

PACS has been developed by a consortium of 
institutes led by MPE (Germany) and including 
UVIE (Austria); KU Leuven, CSL, IMEC (Bel- 
gium); CEA, LAM (France); MPIA (Germany); INAF- 
IFSI/OAA/OAP/OAT, LENS, SISSA (Italy); IAC 
(Spain). This development has been supported by 
the funding agencies BMVIT (Austria), ESA-PRODEX 
(Belgium), CEA/CNES (France), DLR (Germany), 
ASI/INAF (Italy), and CICYT/MCYT (Spain). HIPE 
is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Her- 
schel Science Center, and the HIFI, PACS and SPIRE 
consortia. This work is based [in part] on observations 
made with Herschel, a European Space Agency Corner- 
stone Mission with significant participation by NASA. 
Support for this work was provided by NASA through an 
award issued by JPL/Caltech. This research has made 
use of the NASA/IPAC Extragalactic Database (NED) 
which is operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 
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TABLE 1 

PACS Spectroscopic Observations of NGC 1097 



Obs. ID R.A. Dec. Mode 1 Raster Raster Duration Cycles Repe- On-Source Line 

(J2000) (J2000) Pattern Step (") (s) titions (s raster" 1 ) 



1342189609 


02 


46 


14 


-30 17 00 


WS 


2x2 


18x47 


675 


1 


2 


320 


CII 


157.7 


1342189608 


02 


46 


12 


-30 17 10 


CN 


2x1 


12x12 


13200 


2 


8 


6144 


or 


63.18 


1342189415 


02 


16 


14 


-30 15 04 


WS 


2x2 


12x12 


1644 


1 


1 


160 


CII 


157.7 






// 




ll 


ii 


ll 


II 


ll 


ll 


3 


480 


Nil 


121.9 


1342189412 


02 


46 


23 


-30 17 50 


WS 


2x2 


12x12 


758 


1 


1 


160 


CII 


157.7 


1342189411 


02 


46 


19 


-30 16 28 


WS 


2x2 


12x12 


2262 


1 


1 


160 


CII 


157.7 






// 




ll 


II 


ll 


II 


ll 


ll 


1 


160 


Nil 


121.9 


1342189410 


02 


46 


19 


-30 16 28 


CN 
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12x12 


13833 


1 


1 


704 


CII 


157.7 






// 




ll 


n 


II 


ll 


ll 


II 


1 


704 


OI 


63.18 
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ll 
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ll 


II 


ll 


II 
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Nil 
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Nil 
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1 


9 
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WS 


2x2 


12x47 


945 


1 


3 
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CN 


2x1 


12x12 


7459 
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2 
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Note. — Summary of PACS spectral observations of NGC 1097. Units of right ascension are hours, minutes, and seconds, and units 
of declination are degrees, arcminutes, and arcseconds. When multiple lines were grouped into an observation each line is listed on a 
separate line; only observational parameters that are line specific are given in additional lines of an observation. 

1 CN = Chop-Nod; WS = Wavelength Switching 

2 Parts of these observations were affected by grating malfunctions. Affected raster positions have not been included in combined maps. 



TABLE 2 

PACS Spectroscopic Observations of NGC 4559 



Obs. ID 


R.A. 


Dec. 


Mode 1 


Raster 


Raster 


Duration 


Cycles 


Repe- 


On-Source 


Line 




(J2000) 


(J2000) 




Pattern 


Step (") 


(s) 




titions 


(s raster -1 ) 






1342187229 


12 36 02 


+27 56 26.9 
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2x1 


16x14.5 


6674 


1 


8 


1536 


[O I] 


63.18 


1342187230 
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// 
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[O I] 
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1 
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1 
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ll 


ll 


II 


II 


II 


n 


II 


2 


80 


C II 


157.7 


1342189068 


12 35 52 
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Note. — Summary of PACS spectral observations 


of NGC 4559. Units of 


right ascension arc 


hours, minutes, and seconds, and units 



of declination arc degrees, arcminutes, and arcseconds. When multiple lines were grouped into an observation each line is listed on a 
separate line; only observational parameters that are line specific are given in additional lines of an observation. 

1 CN = Chop-Nod; WS = Wavelength Switching 

2 Parts of these observations were affected by grating malfunctions. Affected raster positions have not been included in combined maps. 



